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ABSTRACT 

We report the results of our multicolor observations of PG 1 1 15+080 with the 1 .5-m telescope 
of the Maidanak Observatory (Uzbekistan, Central Asia) in 2001-2006. Monitoring data in fil- 
ter R spanning the 2004, 2005 and 2006 seasons (76 data points) demonstrate distinct bright- 
ness variations of the source quasar with the total amplitude of almost 0.4 mag. Our R light 
curves have shown image C leading B by 16.4d and image (A1+A2) by 12d that is inconsis- 
tent with the previous estimates obtained by Schechter et al. in 1997 - 24. 7d between B and 
C and 9.4d between (A1-HA2) and C. The new values of time delays in PG 1115+080 must 
result in larger values for the Hubble constant, thus reducing difference between its estimates 
taken from the gravitational lenses and with other methods. Also, we analyzed variability of 
the A2/A1 flux ratio, as well as color changes in the archetypal "fold" lens PG 1115+080. 
We found the A1/A2 flux ratio to grow during 2001-2006 and to be larger at longer wave- 
lengths. In particular, the A2/A1 flux ratio reached 0.85 in filter / in 2006. We also present 
evidence that both the Al and A2 images might have undergone microlensing during 2001- 
2006, with the descending phase for Al and initial phase for A2. We find that the A2/A1 flux 
ratio anomaly in PG 1115 can be well explained both by microlensing and by finite distance 
of the source quasar from the caustic fold. 

Key words: cosmology: gravitational lensing - galaxies: quasars: individual: PG 1 1 15+080. 



1 INTRODUCTION 

Gravitationally lensed quasars are known to potentially provide es- 
timates of the Hubble constant Hq from measurements of the time 
delays between the quasar intrinsic brightness variations seen in 
different quasar images (Refsdal 1964). Since a phenomenon of 
gravitational lensing is controlled by the surface density of the total 
matter (dark plus luminous), it provides a unique possibility both to 
determine the value of Hq and to probe the dark matter content in 
lensing galaxies and along the light paths in the medium between 
the quasar and observer 

By now the time delays have been measured in about 20 grav- 
itationally lensed quasars resulting in the values of Hq that are gen- 
erally noticeably less than the most recent estimate of Ha obtained 
in the HST Hubble Constant Key Project with the use of Cepheids ~ 
Hq = 72+8 km s"' Mpc"' (Freedman et al. 2001). This discrepancy 
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is large enough and, if the Hubble constant is really a universal con- 
stant, needs to be explained. A detailed analysis of the problem of 
divergent Ho estimates inherent in the time delay method, and the 
ways to solve it can be found, e.g., in Keeton & Kochanek (1997), 
Saha & Williams (1997), Kochanek (2002), Kochanek &. Schechter 
(2004) and Schechter (2005), Read et al. (2007), and in many other 
works. 

The main sources of uncertainties in determining Hq are: 

• low accuracy of the time delay estimates caused by poorly 
sampled and insufficiently accurate light curves of quasar compo- 
nents, as well as by microlensing events and, as a rule, by low am- 
plitudes of the quasar intrinsic variability; 

• difference in the values of cosmological constants adopted in 
deriving Hq, 

• invalid models of mass distribution in lensing galaxies. 

The way to reduce the effect of the first source of errors is 
clear enough: more accurate and better sampled light curves of a 
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sufficient duration are needed. A chioice of the cosmological model 
is usually just indicated - this is mostly a question of agreement. 
As to the third item, here the problem of estimating the Hubble 
constant encounters the problem of the dark matter abundance in 
lensing galaxies. 

The problem of determining the Hubble constant from the 
time delay lenses is known to suffer from the so-called central con- 
centration degeneracy, which means that, given the measured time 
delay values, the estimates of the Hubble constant turn out to be 
strongly model-dependent. In particular, models with more cen- 
trally concentrated mass distribution (lower dark matter content) 
provide higher values of H^, more consistent with the results of 
the local Hq measurements than those with lower mass concentra- 
tion towards the center (more dark matter). Moreover, it has long 
been noticed that the time delays are sensitive not only to the total 
radial mass profiles of lensing galaxies, but also to the small pertur- 
bations in the lensing potential (e.g., Blandford & Narayan 1986, 
Witt et al. 2000, Oguri 2007). It is interesting to note that this effect 
has been recently proposed as a new approach to detect dark matter 
substructures in lensing galaxies (Keeton & Moustakas 2009). 

The Hubble constant - central concentration degeneracy is a 
part of the well known total problem of lensing degeneracies: all 
the lensing observables, even if they were determined with zero er- 
rors, are consistent with a variety of the mass distribution laws in 
lensing galaxies. A strategy for solving this non-uniqueness prob- 
lem could be a search through a family of lens models that are 
capable of reproducing the lensing observables (Williams & Saha 
2000, Oguri et al. 2007). Then many models can be run in order to 
infer a probability density for a parameter under investigation, e.g. 
for //fl (Williams & Saha 2000). The most recent studies (Saha et 
al. 2006, Read et al. 2007) have shown that, in such an approach, 
the discrepancy between the Ho value determined from lensing and 
with other methods can be substantially reduced if non parametric 
models for mass reconstruction are used, which can provide much 
broader range of models as compared to the parametric ones. 

In defining priors on the allowed space of lens models, it is 
naturally to assume that lensing galaxies in the time delay lenses are 
similar in their mass profiles to other early-type ellipticals, that are 
presently believed to be close to isothermal and admit the presence 
of the cold dark matter haloes. The isothermal models are also con- 
sistent with stellar dynamics, as well as with the effects of strong 
and weak lensing. 

The quadruply lensed quasars are known to be more promis- 
ing for solving these problems as compared to the two-image lenses 
since they provide more observational constraints to fit the lens 
model. Ten astrometric constraints can be presently regarded as 
measured accurately enough for most systems. This especially con- 
cerns the relative coordinates of quasar images. As to the lensing 
galaxiy, its less accurate coordinates are often the only reliable in- 
formation about the lensing object known from observations, with 
other important characteristics being derived indirectly. This situ- 
ation is inherent, e.g. in PG 1115-1-080 with its faint, 0.31- red- 
shift galaxy. Of other observational constraints, the time delays and 
their ratios are very important. In quadruple lenses, the time delay 
between one of the image pairs is usually used to determine H^, 
while the other ones form the //o-independent time delay ratios to 
constrain the lens model, (Keeton & Kochanek 1997). 

It has long been known that the observed positions of multiple 
quasar macroimages are well predicted by smooth regular models 
of mass distribution in lensing galaxies, while their brightness ra- 
tios are reproduced by such models poorly, (e.g., Kent and Falco 
1988, Kochanek 1991, Mao & Schneider 1998). The first system- 
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Figure 1. PGl 115+080 from observations in filter R with the 1.5-m tele- 
scope of the Maidanak Observatory. The image was obtained by averaging 
of six frames from a series obtained in Febniary 24, 2004, with a subsequent 
Richardson-Lucy processing. 

atic analysis of this problem called "flux ratio anomalies" was made 
by Mao & Schneider (1998), who assumed that the anomalies of 
mutual fluxes of the components in some lenses can be explained 
by the presence of small-scale structures (substructures) in lensing 
galaxies or somewhere near the line of sight. 

A popular model of forming hierarchical structures in the Uni- 
verse with a dominant content of dark matter is currently known to 
poorly explain the observed distribution of matter at small scales. 
In particular, the expected number of satellite galaxies with masses 
of the order of Mq » 10^ Mq remained after the process of hierar- 
chical formation is completed, is an order of magnitude larger than 
a number of dwarf galaxies with such masses actually observed 
within the Local Group (Klypin et al. 1999, Moore et al. 2001). 
One of the solutions of this contradiction is a suggestion that some 
substructures, especially those with low masses, are not luminous. 

Metcalf and Madau (2001) were the first to note that the dark 
matter paradigm can naturally explain existence of substructures 
in galaxies lensing the remote quasars, as proposed by Mao and 
Schneider (1998) to interpret the anomalies of mutual fluxes of 
quasar macroimages, and vice versa, confirmation of substructures 
with masses from 10* Mq to 10^ Mq is capable of removing the con- 
tradiction between the predicted number of the low-mass satellite 
galaxies and that one actually observed. The idea turned out to 
be intriguing and was immediately taken up, (Bradac et al. 2002, 
Chiba 2002, Dalai & Kochanek 2002, Keeton 2001, Metcalf & 
Zhao 2002). Investigation of flux ratio anomalies in gravitationally 
lensed quasars is presently believed to be a powerful tool in solv- 
ing the problem of the dark matter abundance in the Universe. It is 
intensively discussed in numerous recent publications (Congdon & 
Keeton 2005; Keeton et al. 2003, 2005; Kochanek & Dalai 2004; 
Mao et al. 2004; Miranda & Jetzer 2007; Pooley et al. 2006, 2007, 
2009; Morgan et al. 2008). 

In Sec. 3 we report our measurements of the A2/A1 flux ratios 
in filters V, R and / from our data obtained in 2001, 2002 and 2004- 
2006 at the Maidanak Observatory and analyse their behaviors in 
time and in wavelength. In Sec. 4, we analyze the new estimates of 
the time delays between the PG 1115-)-080 images, obtained from 
our monitoring in the R filter during 2004-2006 and reported in 
Vakulik et al. (2009). The new values differ noticeably from those 
reported by Schecter et al. (1997) and Barkana (1997). In Sec. 5, 
we discuss our results and their possible effect on selecting an ade- 
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Figure 2. The light curves of PG 1115+080 Al, A2, B, C from observations in filter R with the 1.5-m telescope of the Maidanak Observatoiy in 2001, 2002, 
2004, 2005 and 2006. 



quale lens model for PG 1 1 15+080 and estimating the value of the 
Hubble constant. 



2 OBSERVATIONS 

The quadruply imaged quasar PG 1115+080 is one of the most 
promising candidates both to investigate the dark matter problem 
and to determine the Ho value from measurements of the time de- 
lays between the image components. The source with a redshift 
of zs = 1.722 is lensed by a galaxy with za = 0.31 (Henry & 
Heasley 1986, Christian et al. 1987, Tonry 1998), which forms four 
quasar images, with an image pair A 1 and A2 bracketing the critical 
curve very close to each other. It is the second gravitationally lensed 
quasar discovered over a quarter of century ago, at first as a triple 
quasar (Weymann et al. 1980). Hege et al. (1980) were the first to 
resolve the brightest image component into two images separated 
by 0.48 arcsec. Further observations (Young et al. 1981, Vanderri- 
est et al. 1986, Christian et al. 1987, Ki-istian et al. 1993, Courbin 
et al. 1997) have provided positions of quasar images and infor- 
mation about the lensing object, which allowed construction of a 
macrolens model (e.g., Keeton, Kochanek & Seljak 1997). In par- 
ticular, Keeton & Kochanek (1997) have shown that the observed 
quasar image positions and fluxes and the galaxy position can be 
fit well by an ellipsoidal galaxy with an external shear rather than 
by only an ellipsoidal galaxy, or by a circular galaxy with an exter- 
nal shear. They noted that a group of nearby galaxies detected by 
Young et al. (1981) could provide the needed external shear. 

Observations of PG 1115+080 were started at the 1.5-meter 
telescope of the high-altitude Maidanak Observatory (Central Asia, 
Uzbekistan) in 2001. An image scale of 0.26"/pix was available at 
the f/8 focal plane with a scientific BroCam CCD camera having 
a SITe ST 005A 2030 x 800 chip. The CCD images were usually 
taken in series consisting of 2 to 10 frames for the R filter and of 
2 to 6 frames for V and /. To provide higher photometric accu- 
racy, we averaged the values of magnitudes estimated from indi- 



vidual frames. The seeing varied from 0.75" to 1.3" [the full width 
at half-maximum (FWHM) of images of the reference stars B and C 
according to designation by Vanderriest et al. (1986)]. The analysis 
of photometry shows no significant dependence of the photometry 
errors on seeing, excepting the FWHM noticably exceeding 1.3". 
Occasional frames with such values of the FWHM were excluded 
from processing. 

In Fig. 1 we show one of the best images of PG 1115+080 
obtained through the R filter. For better view, the image was re- 
stored with an algorithm similar to that proposed by Richardson 
(1972) in optics and independently by Lucy (1974) in astronomy 
(the Richardson-Lucy iterative method). 

Our algorithm for photometric image processing is similar to 
that applied to Q2237+0305 and described in detail by Vakulik et 
al. (2004). The light curves of PG 1115+080A1,A2,B,C in filter 
for the time period from April 2001 to lune 2006 are shown in Fig. 
2. The photometry in the V, R and / filters is presented in Tables 
6-8. 

Unfortunately, observations were not carried out in 2003, and 
the data are very scanty for the 2001 and 2002 seasons in all the 
three filters. The most numerous data were obtained in filter R, 
especially in 2004 (23 nights), 2005 (27 nights) and in 2006 (24 
nights). The data demonstrate noticeable variations of the quasar 
brightness, with the total amplitude reaching approximately 0.4 
mag in 2004-2006, and smaller amplitudes of about 0.05 mag on 
a time-scale of two months, which are clearly seen in all the four 
light curves in 2004. The time delays between the light curves of 
the C and B, C and Al (or A2) images can be easily seen from a 
simple visual inspection of the R light curves, therefore, the data 
obtained in 2004-2006 in filter R seem to have good prospects for 
obtaining reliable estimates of the time delays in PG 1 1 15+080. 

Thanks to the spatially resolved photometry of the Al and A2 
image pair in filters V, R and /, our data have made it possible to 
measure flux ratios for these components for five seasons of ob- 
servations, and to study their behavior in time and in wavelength. 
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Table 1. Estimates of the A2/A1 brightness ratios in PG 1 1 15+030 for the time period 1980-2006 from all available data. 



Date 


A2/A1 
flux ratio 


Spectral 
range 


Instrument 


Reference 


1980 June 


0.83 


V 


MMT 


Hege et al. (1980) 


1981 April 30 


1.0 


B 


CFHT 


Vanderriest(1986) 


1983 March 8 


1.0 


V 


_ " _ 


_ " _ 


1984 March 26 


0.95 


B 




_ " _ 


1985 March 16 


0.75 


B 




_ " _ 


1985 March 19 


0.79 


V 




_ " _ 


1986 February 19 


0.79 


V 


CFHT 


Christian (1987) 




0.8 


R 


_ " _ 


_ " _ 




0.79 


B 


_ " _ 


_ " _ 


1989 April 


0.68 


I 


CFHT 


Schechter(1993) 


1991 March 3 


0.66 


V 


HST 


Kristian (1993) 


_ " _ 


0.7 


I 


-"- 


_ " _ 


1992 April 


0.67 


I 


Hiltner 


Schechter 1993 


-"- 


0.69 


V 


-"- 


_ " _ 


_ " _ 


0.72 


I 


CTIO 


_ " _ 


_ " _ 


0.68 


V 


-"- 


_ " _ 


1993 April 


0.69 


I 


Hiltner 


_ " _ 


_ " _ 


0.63 


V 


-"- 


_ " _ 


1995 December 20 


0.66 


V 


Magellan 


Pooley et al. (2006) 


1996 June 7 


0.68 


I 


NOT 


Courbin et al. (1997) 


1997 November 17 


0.64 


H 


HST 


Impey et al. (1998) 


? 


0.52 


V 


HST 


Morgan et al. (2008) 


? 


0.67 


I 


" 


_ " _ 


? 


0.63 


H 




_ " _ 


2001 March 26 


0.66 


V 


Magellan 


Pooley et al. (2006) 


2001 April 20-27 


0.74 


V 


1.5m(Maidanak) 


This work 




0.67 


R 


-"- 


_ " _ 




0.72 


I 


_ " _ 


_ " _ 


2002 March 


0.76 


V 




_ " _ 


_ " _ 


0.71 


R 




_ " _ 




0.72 


I 






2004 February 22 


0.81 


Sloan i' 


Magellan 


Pooley et al.(2006) 


2004 May 5-6 


0.93 


ll.67fim 


Subaru 


Chiba et al. (2005) 


2004 Jan. 17-June 8 


0.79 


V 


1.5m Maidanak 


This work 




0.81 


R 






2004 Apr. 11 -June 8 


0.83 


I 






2005 June 07 


0.81 


Sloan i' 


Magellan 


Pooley et al. (2006) 


2006 Jan.5-Apr.l5 


0.8 


V 


1.5m Maidanak 


This work 


2006 Jan.5-June 2 


0.83 


R 






2006 Jan.5-Apr.l5 


0.85 


I 







As is noted above, deviations of flux ratios in quasar macroimages 
from the theoretical predictions (flux ratio anomalies) are presently 
believed to be diagnostic for detection of substructures in lensing 
galaxies, which may represent the dark matter. 



3 THE A2/A1 FLUX RATIOS 

The idea to detect substructures in lensing galaxies using the 
anomalies of flux ratios is based on fundamental relationships be- 
tween coordinates and magnifications of the quasar images, which 
result from the general lens equation (Schneider, Ehler & Falco, 
1992). These relationships have been obtained for the first time by 
Schneider and Weiss (1992) and Mao (1992) for several "smooth" 
distributions of lensing potential. In principle, the lens equation is 
capable of providing six independent relationships between the co- 
ordinates and magnifications for a quadruple lens, but only one of 
them can be checked with the data of observations. This is the well- 
known magnification sum rule for a source within a macrocaustic 



cusp, when three close images emerge: magnification of the central 
image must be equal to the sum of magnifications of two outer im- 
ages (Schneider & Weiss 1992). When the source lies near a caus- 
tic fold, two images of the same brightness must arise (Keeton, 
Gaudi and Fetters 2005, KGP hereafter). Since the absolute val- 
ues of magnifications in macroimages are unknown (the unlensed 
quasar cannot be observed), Mao & Schneider (1998) proposed to 
use the dimensionless quantities 

^ ^ ly"il - ly"2l + lj"3l ^ Fi - F2 + F3 
I//1I + I//2I + 1^/31 " F1+F2 + F3 

for three images emerging when the source is in a caustic cusp, and 

\li,„\ - 1^,1 F,„ - F,, 
Kfoid - ; — — — : - — — — rr 

\fJm\ + F„, + F., 

for the case when the source is at the caustic fold. Indices m and 
i- in the second expression denote the images at the minimum and 
saddle points of the Fermat surface. 
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Figure 3. A history of the A2/A1 flux ratios in PG 1 1 15+080 from the data 
in filters V and / as listed in Table 1 . 

Ideally, cusp relation Rcusp = and fold relation Rfoid = 
hold only when the source lies exactly at the caustic cusp or 
fold, respectively. In real lenses these relations hold only approx- 
imately. KGP(2003, 2005) fulfilled a detailed study of asymptotic 
behaviours of the cusp and fold relations and calculated probability 
distributions of Rfou values for several smooth lens models. The 
value of deviation of Rcusp and Rfdd from zero can be regarded as 
a measure of probability for the lensing potential to have substru- 
tures on scales smaller than the separation between the closest im- 
ages (KGP 2003, 2005). KGP (2005) warn, however, that for fold 
lenses, the observed violation of the fold relation may just mean 
that the source is far enough from a caustic fold. 

It should be noted that, in principle, the observed anomalies 
of brightness ratios in images of gravitationally lensed quasars can 
be explained by other factors, such as microlensing by compact 
bodies and the effects of propagation phenomena in the interstel- 
lar medium (extinction and scattering, scintillations). These factors 
are studied in details by Kochanek & Dalai (2004). They concluded 
that substructures of cold dark matter is the best explanation for the 
flux ratio anomalies in some quadruply lensed quasars. They re- 
minded also that, as was stated for the first time by Mao & Schnei- 
der (1998), the fluxes of highly magnified saddle images are very 
sensitive to small gravitational perturbations as compared to low- 
magnification images and, even more importantly, these perturba- 
tions bias the fluxes towards demagnification, as was also noted by 
Schechter & Wambsganss (2002). 



3.1 Behavior of the flux ratios in time 

The A1+A2 image pair in PG 1 1 15+080 consists of a highly mag- 
nified minimum point image (Al) and saddle point image (A2) sit- 
uated symmetrically with respect to the fold caustic very close to 
each other. According to theoretical expectations (e.g. Schneider et 
al. 1992), the ratio of their fluxes must be close to 1. 

There are numerous measurements of the A2/A1 brightness 
ratio in PG 1 1 15+080 made at different spectral ranges and at dif- 
ferent epochs since 1980, which we tried to assemble in Table 1. 
Some of these data have been used by Pooley et al. (2009) to anal- 
yse a long-term history of the A2/A1 variations in the optical band 
and to compare with the X-ray data, (see, e.g., Pooley et al. 2006, 
2007, 2009). Fig. 2 from the work by Pooley et al. (2009) demon- 
strates changes in the A2/A1 optical flux ratio during the time pe- 
riod from 1980 to 2008, and much more dramatic changes of this 
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ratio in X-rays. Pooley et al. (2009) noted that, according to all 
the observations since the system discovery, the A2/A1 flux ratio 
varied within 0.65-0.85. They did not specify the optical spectral 
bands for the data in their Fig. 2, however, but it looks like they are 
for filter V. 

In our Table 1, the telescopes and filter bands are indicated for 
all estimates of the A2/A1 flux ratios. The table does not contain the 
results of the Chandra X-rays observations, which exhibited strong 
flux ratio anomaly and can be found, e.g., in Pooley et al. (2006, 
2007, 2008). The data in Table 1 are also supplemented by the esti- 
mates of the A2/A1 flux ratios obtained from our photometry (Ta- 
bles 5-7). Flux ratios in filters V and / from Table 1 are displayed 
in Fig. 3. These flux ratios behave similarly in time for both filters, 
and in general features resemble fig. 2 from Pooley et al. (2008). 
Based upon their fig. 2, Pooley et al. (2008) argue that the optical 
flux ratio anomaly in PG 1 1 15+080 is slight and "nearly constant in 
time". Our analysis described below have shown, however, that it is 
not quite so. Our Fig. 3, where the available previous data in V and 
/ are supplemented by our measurements, shows that variations of 
the A2/A1 flux ratio in time are indeed rather small and slow. How- 
ever, even if we exclude a marginal value for the date 1983, March 8 
(Vanderriest et al. 1986), which equals 1 with the uncertainty of 0. 1, 
we will have the A2/A1 flux ratio varying in some regular manner 
with the amplitude of about 0.15 during the last 25 years. Some- 
where between 1991 and 1996, the ratio reached its minimal value 
of about 0.65 in filter V, and increased up to 0.8 by 2006. It should 
be noted that the fact that A2/A1 flux ratio varies in time is in it- 
self an argument in favour of microlensing as the main reason for 
the anomalous flux ratio in PGl 115+080. Also, it should be men- 
tioned that the A2/A1 flux ratio is slightly but steadily higher in 
filter / than in V. 

To determine which component (or components) exactly un- 
derwent microlensing, we addressed only our data as more homo- 
geneous ones, and analysed behaviors in time of the long-term con- 
stituents of the A2-A1, C-Al, B-C and C-A2 magnitude differences 
for filters R and /. These difference light curves in filter R are shown 
in Fig. 4. We did not correct the individual light curves for the time 
delays, which are small as compared to the characteristic time-scale 
of quasar flux variations. This might result only in some increase 
of the data points scatter with respect to the approximating curves, 
which are the second-order polynomials in Fig.4. 

The largest decrease of the magnitude difference is for the 
A2-A1 image pair - about 0.23 mag during 2001-2005. Pair C- 
Al shows an almost linear decrease of the magnitude difference 
in time, with only 0.12 mag during 2001-2005. Since the mutual 
brightness of images B and C was almost invariable in 2001-2006, 
one might conclude that it is an image Al that became fainter 
during this time period. But the C-A2 magnitude difference curve 
shows however, that, in addition to the obvious dimming of image 
Al, brightening of image A2 makes a certain contribution to the 
decrease of the A2-A1 magnitude difference in 2001-2005. 

Therefore, we can conclude that a decay of Al and brighten- 
ing of A2 took place simultaneously in PG 1 1 15+080 during 2001- 
2006. We may also conclude that it is the Al image that under- 
went microlensing in the previous years, with the maximum near 
1992-1995, as seen from Fig. 2, and the final phase in 2006 or, 
perhaps, later. With the previous data taken into account (see Fig. 
3 and Table 1), the total time-scale of the 0.3-magnitude event is 
about 25 years. Image A2 underwent microlensing as well, with 
its rising branch occurring in 2001-2005. The brightening of image 
A2 reached about 0. 14 mag during this time period, while the to- 
tal brightening in the whole event may be larger. In calculation of 
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Figure 4. Behaviors of the C-Al, C-A2, B-C and A2-A1 magnitude differences in time from the results of our photometry in filter R\ approximation by the 
second-order polynomials is shown. 



Table 2. Flux ratios in PG 1115+080 as predicted by the most recent lens 
models and determined from the results of our photometry in 2006 (filter 
I); the uncertainty of our flux ratio estimates is 0.02 for all ratios. 



Lens model 


A2/A1 


B/Al 


C/Al 


B/C 


Chiba (2002) 


0.92 


0.22 


0.28 


0.8 


Chiba (2005) 


0.92 


0.22 


0.28 


0.79 


Pooley(2006) 


0.96 




0.26 


0.67 


Pooley(2007) 


0.92 


0.21 


0.27 


0.78 


This work 


0.85 


0.19 


0.29 


0.68 



the time delays, more subtle variations of the magnitude differences 
during 2004-2006 were found, (see Sec. 4.1). 

It should be noted that our results are well consistent with 
measurements of the A1-A2 magnitude difference presented in 
Morgan et al. (2008), who reported approximately 0.2-mag growth 
of this quantity during 2001-2006. However, they do not present 
the magnitude differences between other images and Al or A2 sep- 
arately, which has led us to a conclusion about the final phase of 
microlensing in image Al and, seemingly, the initial phase of a 
microlensing event in image A2. This conclusion is also indirectly 
confirmed by the results of Pooley et al. (2009), who reported a dra- 
matic rise in the X-ray flux from image A2 between 2001 and 2008. 
Larger microlensing amplitudes at shorter wavelengths are often 
detected for many lensed quasars and are known to be naturally ex- 
plained by smaller effective sizes of quasars at shorter wavelengths. 

The observed time-scales and amplitudes of the microlensing 
brightness fluctuations are known to depend on the relative veloc- 
ity of a quasar and lensing galaxy, and on the relationship between 
the source size and the Einstein ring radius of a microlens. For 
PGl 115+080, the expected duration of a microlensing event is es- 



timated to be of the order of 10 to 20 years for the subsolar mass 
microlens (Chiba 2005), well consistent with that in image Al ob- 
served in 1980-2006. 

Thus, if our interpretation of the observed brightness varia- 
tions in the Al, A2, B and C images is valid, then the A2/A1 flux 
ratio would be expected to approach its undisturbed value in a few 
years, unless a new event takes place in at least one image. Our es- 
timates of the A2/A1, B/Al, C/Al and B/C flux ratios calculated 
from the photometry data of 2006, are presented in Table 2 along 
with the model predictions by Chiba (2002), Chiba et al. (2005), 
Pooley et al. (2006, 2007). The ratios for filter / are presented, 
where the effect of microlensing is expected to be minimal as com- 
pared to V and R. As is seen from Table 2, the A1/A2 flux ratio 
is still less than predicted by the most recent lens models. How- 
ever, when expressed in terms of Rf„u, it would equal 0.08, which 
means that, according to simulations of KGP (2005), this flux ratio 
is within a region admissible by a smooth lensing potential model 
for the finite source distances from the caustic fold, that is, it is not 
anomalous in the sense implied by Mao&Schneider (1998). 

3.2 Color Changes in PG1115-H080 

We have also made use of our multicolor observations to analyze 
behaviors of the V - I color indices of image components, which 
were shown to be indicative of the microlensing nature of bright- 
ness changes for at least the Q2237+0305 quasar (Vakulik et al. 
2004). Since only 16 data points were available to build the V - / 
versus R dependency for each image component, we did not build 
them for the components separately, but combined the data into two 
sets, for A1+A2 and B+C image pairs. The resulting diagrams are 
presented in Fig. 5. In order to eliminate the magnitude difference 
and possible permanent color difference between the components 
in each pair, we shifted the data points in both plots along the R 
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Figure 5. V - I versus R diagrams for image A (upper panel) and B,C 
(bottom), illustrating a statistical dependence between the color indices and 
magnitudes. The regression line slopes differ insignificantly for the A1,A2 
and B,C image pairs and equal 0.29 and 0.28, with the congelation indices 
0.61 and 0.56, respectively. 
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Figure 6. A history of the long-term variations of color indices V - / in 
PGl 1 15+080. Each point is a result of averaging ofV — I values within the 
corresponding season; the midpoints of every season are along the horizon- 
tal axis. 



axis by the values of the mean magnitude differences between the 
components. The values of V - I grew with the growth of the R 
magnitude in both diagrams, with the regression line slopes of 0.28- 
0.29. This is qualitatively consistent with the total 0.4-mag fading 
of the PGll 15-1-080 quasar during 2001-2005: according to numer- 
ous observations, there is a common tendency for many quasars to 
become bluer at their bright phases (see, e.g. Giveon et al. 1999 and 
references therein). In particular, Giveon et al. (1999) presented the 
A(B - R) versus AS and A(S - R) versus AR diagrams built for a 
subset of 21 quasars from their Palomar Green sample consisted of 
42 quasars. Their diagrams show a significant correlation between 
the color and magnitude variations, with the regression line slopes 
of 0.25-0.27. 

The regression line slopes for the diagrams in our Fig. 5 are 
also close to that reported for Q2237-I-0305 in one of our previous 
publications (Vakulik et al. 2004). There is an important difference 
between the two quadruple lenses, however: the Q2237+0305 light 
curves are strongly dominated by microlensing events, while in 
PGl 1 15-1-080, a contribution from microlensing activity is small as 
compared to the quasar intrinsic variability. We do not see any sig- 
nificant difference between the two diagrams in Fig. 5, though the 
contributions from microlensing for these image pairs are different. 
Therefore, the diagrams in Fig. 5 should be referred to characteris- 
tics of the PGl 1 15-1-080 quasar variability rather than to microlens- 
ing variability, in contrast to Q2237, where they result mostly from 
microlensing. 

Also, we tried to analyze the long-term history of color 
changes in PGl 1 15-1-080 and plotted the values of V - I averaged 
within every season as functions of the corresponding time mo- 
ments, which are the midpoints of the seasons (Fig. 6). It should be 
remembered that the values of color indices in this plot are accurate 
to only ~ ±0.025mag, therefore. Fig. 6 illustrates their behaviors in 
time qualitatively rather than quantitatively. Nevertheless, though 
the data points in this plot are rather scattered and often overlap 
giving a rather intricate pattern, the general features in behaviors 
of color indices are evident. To make them clearer, we showed the 
corresponding linear approximations in Fig. 6. First we notice that 
the V - I color indices increase in time for all the four image com- 
ponents. This conclusion seems to be valid, while the differences in 



gradients of the V - I colors with time are hardly veridical. In gen- 
eral. Fig. 6 confirms the dominating contribution from the quasar in- 
trinsic variability over microlensing activity in PG 1 1 15+080 light 
curves, as was indicated in the section above. Another significant 
conclusion from Fig. 6, is that the V - I color index of image B re- 
mains permanently smaller that those of the three other ones, while 
demonstrating similar increase as well. We fail to explain why an 
image that lies the smallest distance from the lensing galaxy turned 
out to have the bluest color index. 



4 ANALYSIS OF THE NEW ESTIMATES OF TIME 
DELAYS IN PG 1115-H080 

The time delays in PG 1 1 15+080 were determined for the first time 
by Schechter et al. (1997) to be 23.7 ± 3.4 days between B and 
C, and 9.4 ± 3.4 days between A1+A2 and C (image C is lead- 
ing). Barkana (1997) re-analyzed their data using another algorithm 
and reported 25^3 g days for the time delay between B and C, and 
this is quite consistent with 23.7 ± 3.4 days from Schechter et al. 
(1997). But the other time delays, and hence the time delay ra- 
tio r^BC = TacItba differ significantly: rj^BC — 1-13;^q |7 as calcu- 
lated by Barkana (1997) and 0.7 ± 0.3 according to Schechter et 
al. (1997). Since 1997, just these values, either the first or the sec- 
ond ones, were being used to constrain the PGl 1 15+080 model and 
determine the Hubble constant. 

Determination of the time delays has generated a flow of mod- 
els for the system, (Schechter et al. 1997, Keeton & Kochanek 
1997, Courbin et al. 1997, Impey et al. 1998, Saha & Williams 
2001, Kochanek, Keeton & McLeod 2001, Zhao & Pronk 2001, 
Chiba 2002, Treu & Koopmans 2002, Yoo et al. 2005, 2006, Poo- 
ley et al. 2006, Miranda & Jetzer 2007), all illustrating how strongly 
the estimated value of Hq depends on the adopted mass profiles of 
the lens galaxy for the given values of time delays. 

The detailed analysis of the uncertainties in determining 
Hubble constant from the time delay lenses can be found in, 
e.g., Kochanek & Schechter (2004), Schechter (2005), Kochanek 
(2002), where the paths to eliminate or at least to lessen the un- 
certainties have been also outlined. Kochanek & Schechter (2004) 
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Table 3. The time delays (days) for PGl 115+080 as reported earlier 
and in our previous work, (Vakulik et al. 2009). 



Author 




T"AC 


tbc 


Schechter (1997) 


14.3 ± 3.4 


9.4 ±3.4 


23.7 ± 3.4 


Barkana (1997) 


1 1 7+-' 


13.3!2-0 


25.0!|3 


Vakulik (2009) 


4 4+3.2 


12.0!|§ 


16.4!^J 



indicated, in particular, the importance of improving the accuracy 
of time delays for PG 1 1 15+080. 

The R hght curves taken from the 2004-2006 data (Fig. 2) 
clearly demonstrate their applicability to determine the time de- 
lays. As compared to the data used by Schechter et al. (1997) and 
Barkana (1997), we were lucky to detect the quasar brightness vari- 
ation with an amplitude of almost a factor of three larger, and with 
rather well-sampled data points within every season of observa- 
tions. In addition, the accuracy of our photometry has made it pos- 
sible to confidently detect flux variations with an amplitude as small 
as 0.05 mag that can be seen in the data of 2004. 

The methodology to determine the time delays in pairs is sim- 
ple enough and obvious. A common feature of all known methods 
of time delay measurements is the use, in one way or another, of the 
cross-correlation maximum or mutual dispersion minimum criteria, 
while they may differ in the algorithms of the initial data interpola- 
tion. 

Analysis of the light curves of quasar images in pairs can also 
be applied when a lens consists of more than two images. To deter- 
mine the time delays from the light curves shown in Fig. 2, we used 
another approach however, as described in more detail in our pre- 
vious works (Vakulik et al. 2006, 2009). Here we shall only remind 
the fundamentals of our approach. 

We determined the source light curve from a joint analysis 
of light curves of all image components. The individual time de- 
lays for pairs of images can be then determined with respect to this 
model source light curve jointly from the corresponding system of 
equations. 

Since, according to predictions of all lens models and to mea- 
surements in X-rays by Dai et al. (2001) and Chartas et al. (2004), 
the time delay between images Al and A2 does not exceed a small 
fraction of the day, their fluxes were summed to form a single curve, 
which we call the A light curve. Thus, we may write the following 
functional for three light curves: 

0(Af,To,T|,T2) = — ) ) — i i— '—, 3) 

where mj{tj) are the data points in the light curve of the jth image 
at the time moments f, ; dnij and Tj are the shifts of a corresponding 
light curve in stellar magnitude and in time, respectively, W is a 
number of points in the light curves, Af is the parameter of the 
approximating function /(/,, A?, t^), and (T^(?,) are the photometry 
errors. 

We adopted dmo = and tq = in our calculations, that is, we 
fitted the light curves of B and C to the A light curve, and thus, dnii 
and dm2 are the magnitude differences A-B and A-C, respectively. 
At given values of ti and T2, we minimize <I>(A/, Ti , T2) in dmj and 
in coefficients of the approximating function. The values of mini- 
mum of <I)(Af, Ti,T2) were being looked for at a rectangular mesh 
Ti , T2 with a step of 0.5 d in preliminary calculations, and of 0.2 d 



Table 4. Time delay ratios rAc/'''BA and tacI'^bc for PGl 115+080 as 
predicted by several lens models (the upper part of the table) and deter- 
mined from the existing measurements of the time delays for the system 
(the last three lines). 



Author 


'^ACI'T'BA 




Schechter et al. (1997) 


1.33-1.80 


0.57-0.64 


Keeton&Kochanek (1997) 


1.35-1.47 




Impey et al. (1998) 


1.3 




Chartas et al. (2004) 


1.3 


0.56 


Keeton et al. (2009) 


1.54 


0.61 


Schechter et al. (1997) 


0.66 


0.40 


Barkana (1997) 


1.14 


0.53 


Vakulik et al. (2009) 


2.73 


0.73 



at a final stage. The values of Ti , T2 corresponding to the minimal 
value of "tCA?, ti , T2) were adopted as the estimates of the time de- 
lays Tba and t^^. Fig. 7 shows a distribution of 0(A?, Ti , T2) - <!>„,■„ 
in the space of parameters Tba and Tac calculated for parameter 
A? = 0.12 years. Thus, our estimates of the time delays that can be 
read out at the t^c and Tbc axes against the centre of contours in 
Fig. 7, are Tba = 4.4, Tac = 12.0 days. The time delay Tbc is not 
an independent quantity in our method, and can be determined as a 
linear combination Tbc = Tba + Tac, that is, Tbc = 16.4 days. 

To test our method for robustness and absence of systematics, 
and to estimate the accuracy inherent in our time delay measure- 
ments, we fulfilled a numerical simulation as described in detail in 
(Vakulik et al. 2009). The simulated light curves of the components 
were obtained by shifting the approximating curve /(f,, Af, r^) by 
the proper time delays ri , T2 and magnitude differences, and by 
adding random quantities to imitate the photometry errors. We sim- 
ulated 2000 light curves synthesized as described above, and cal- 
culated the resulting time delays using the procedure, which was 
exactly the same as in the analysis of the actual light curves. The 
results of simulations were used to build the distribution functions 
for errors and to estimate the 95-percent confidence intervals. 

The final values of the time delays and the corresponding un- 
certainties are presented in Table 3, where they can be compared 
with the estimates reported by Schechter et al. (1997) and Barkana 
(1997). 

It is interesting to note that using only the data of 2004, where 
a small-amplitude turn-over in the light curves is detected, we ob- 
tained Tba = 5.0 days, Tac = 9.4 days, and Tbc = 14.4days, consis- 
tent with the estimates obtained from the whole data set. However, 
simulation of errors for only the data of 2004 demonstrates notice- 
ably larger uncertainties, as compared to those calculated from the 
entire light curve. 

The light curves of images A, B and C shifted by the corre- 
sponding time delays and reduced to image A in magnitude are 
shown in figure 2 from our previous paper (Vakulik et al 2009) for 
the approximating function parameter At = 0.12 years. As is seen 
from this picture, the data points for all the three images are very 
well consistent with each other and with the approximating curve. 

Thus we obtained the time delay values, which differ notice- 
ably from those reported by Schechter et al. (1997) and Barkana 
(1997) and used in a variety of models of many authors to de- 
rive the Hubble constant value. The largest differences are for Tbc 
and Tba'- our estimate of Tbc is a factor of 1.5 smaller, while for 
Tba, it is almost three times smaller as compared to the results 
of Schechter (1997) and Barkana (1997). Meanwhile, our values 
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Table 5. Time delays as predicted by the lens models calculated by Schechter et al.(1997) and the values of the Hubble constant Hq obtained from comparison 
of these time delays with those obtained by Schechter et al. (1997), columns 2-4; the Hq values calculated for the same lens models with the time delays 
determined in this work (columns 5-7). 



Lens models, T(days) and //o(Schechter et al. 1997 ) Hq with r^c ™d '''bc from this work 



Model 


TAcidays) 


TBc(days) 


//o(kms-'Mpc-') 


Ho (TAc = 12.0'') 


HoiTBC = 16.4'') 


//o(rnean) 


PMXS 


12.5 


19.9 


84 


104 


121 


113 


ISXS 


6.6 


10.4 


44 


55 


63 


59 


ISEP 


9.7 


15.1 


64 


81 


92 


86 


ISIS 


5.6 


9.7 


41 


47 


59 


53 


ISIS+ 


5.7 


10 


42 


48 


61 


54 




11.2 11.6 12 12.4 12.8 

\c (days) 



Figure 7. Distribution of <I>(A?, ti , T2) - "t;,,,,, in the space of parameters 
TAC and Tab- The innermost contour corresponds to <D(A;, ti,T2) - 
equalling 0.0001, with every next level twice as much than preceding. 

of Tac are rather similar to those of Schechter and, especially, of 
Barkana. 

As is noted in Introduction, the time delay between one of 
the image pairs, say, Tbc, can be used to determine Hq, while the 
time delay ratio r^gc = Tac I^ba is independent of the H^, value and 
can be used to constrain the lens model. Most of the PGl 1 15-H080 
macromodels are consistent in predicting vabc to within 0. 15. As far 
as can be expected from Table 3, the three measurements of time 
delays do not provide the time delay ratios consistent with each 
other and with model predictions. This can be seen from Table 4, 
where we collected several model predictions for vabc together with 
the measurements presented by Schechter et al. (1997), Barkana 
(1997) and Vakulik et al. (2009). Also, we presented here the time 
delay ratios vcba = Tac/tbc, which are connected with rABC by a 
relationship rcsA = ''abc/(1 + '"abc)- These quantities demonstrate 
better agreement between model predictions and measurements. 

The time delay ratios calculated from measurements of 
Schechter and Barkana are seen to be lower as compared to the 
model predictions, while our measurements provide the estimates 
of this quantities exceeding the model predictions, especially for 
^ABC = Tac/tba, which is as large as 2.73 from our data. One should 
admit that such discrepancy is too large, since the largest rABC we 
have found in the literature is that calculated by Schechter et al. 
(1997) for their isothermal ellipsoid model - r^gc = 1-8. 

The reason for this becomes qualitatively clear when address- 
ing the data in Table 3: the shortest time delay Tba is measured with 



almost the same absolute error as the longest one, Tbc, that is, of all 
the three time delays, Tba has the largest relative error. Therefore, 
we are far from arguing our value of Tba to be more trustworthy that 
those obtained by Schechter and Barkana. We regard, however, that 
the values of Tbc and Tac are more reliable and trustworthy. Also, it 
should be noted that none of the macrolens models we could found 
in the literature predicts the values of Tbc larger than 18d, - the 
value 19. 9d from the unrealistic point-mass model in Schechter et 
al. (1997) may hardly be taken into account. 

In the framework of the present publication, we did not intend 
to either propose an extra exotic lens model or recalculate the most 
popular ones to derive the new estimate of the Hubble constant, 
but instead, we have made use of the results of Schechter et al. 
(1997), who calculated five models for the gravitational potential 
of PGl 1 15+080. The time delays tac and Tbc as predicted by their 
five models for D = 1 and Hg = 100km s^'Mpc"', and the corre- 
sponding estimates of Hq obtained with their time delays are shown 
in Table 5 (columns 2-4). We remind that, according to Schechter 
et al. (1997), model PMXS means a point mass with external shear, 
the ISXS model is an isothermal sphere with external shear, and 
ISEP is an isothermal elliptical potential. The ISIS model uses a 
second isothermal sphere for a group of galaxies at approximately 
the same redshift as the main lensing galaxy (Young et al. 1981; 
Henry & Heasley 1986) to represent shear. In the ISIS+ model, the 
uncertainty in the galaxy position is not regarded to be negligible, 
as in the ISXS and ISIS, but its coordinates were taken as two ad- 
ditional free parameters. 

Columns 5 to 7 of Table 5 contain the Hq values estimated 
for the Schechter et al. models with our values of the time delays 
for image pairs AC and BC separately (column 5 and 6), and the 
average between the pairs (column 7). We may conclude that, as 
could be expected, our new estimates of Tbc and Tac provide higher 
values of Hubble constant, which are closer to the most recent value 
obtained in the HST Key Project from observations of Cepheids 
(Freedman et al. 2001). 



5 DISCUSSION AND CONCLUSIONS 

From our new data of excellent quality, we had a possibility to dis- 
criminate between the quasar intrinsic and microlensing brightness 
and color variations for PG 1 1 15+080, and to obtain new estimates 
of the time delays. It may immediately be seen from our R light 
curves, especially for 2004-6 seasons, that significant brightness 
fluctuations with amplitudes exceeding the typical error bars of the 
data points, have been detected. This allows us to recompute the 
values of Hq calculated with the previous estimates of the time de- 
lays for this gravitational lens system. 



© ... RAS, MNRAS 000,[T]^?? 



10 Tsvetkova VS. et al. 



In particular, we find the following: 

• We have studied behaviors of brightness ratios of all the com- 
ponents both in time and in wavelength. We report microlensing 
in Al with an amplitude of about 0.3 mag in filter R on a. time- 
scale of 25 years. The magnification peak in Al took place in 
1992-1995, with the subsequent fading in 2001-2006. The image 
Al flux may be expected to reach its undisturbed value by 2006 
or later. A microlensing event was apparently observed in image 
A2 as well, with its rising branch in 2001-2005, when image A2 
brightened by approximately 0.15 mag. The time scales and am- 
plitudes of both events are consistent with those predicted for this 
object for the Solar mass microlenses (Schechter 2004). We also 
notice that a similar microlensing amplitude was detected for an- 
other quasar with similar redshift distances to the quasar and lens, - 
for the First Lens Q0957-I-561 (Schild & Smith 1991). In fitting the 
2004-2006 data points to the approximating function, very subtle 
signs of microlensing have also been found in image B. 

• Therefore, deviation of the observed A2/A1 flux ratio from 
that predicted by most of the lens models can be well explained 
by microlensing events. An additional contribution to the flux ratio 
"anomaly" may be expected from the source position with respect 
to the caustic fold: when expressed in terms of Rfoid (Eq. 2), the 
brightness difference between Al and A2 would equal 0.08, which 
means that, according to simulations of KGP (2005), it is within 
a region admissible by a smooth lensing potential model, that is, 
it is not anomalous in the sense implied first by Mao & Schneider 
(1998). 

• We have made use of observations in other filters available for 
some dates to analyse behaviors of color indices of the images. The 
V - 1 versus R diagrams built for pairs A1+A2 and B+C demon- 
strate the known tendency of quasars to become bluer at their bright 
phases, with no signs of any contribution from microlensing: the 
diagrams for both image pairs are nearly identical. This can be ex- 
plained either by poor statistics (the data in all the three filters are 
available for only 16 nights, providing rather poor correlation as in- 
dicated in Fig. 5 caption), or by small amplitudes of the microlens- 
ing events under consideration, or by both reasons. 

• An interesting feature of the behaviours of color indices 
should be noted. While all images demonstrate growth of their 
color indices in time, the V -I color indices of image B are slightly 
but steadily less than those of other images. This is rather unex- 
pected, since image B is located the closest distance to the main 
lensing galaxy. 

• The time delays for PGl 1 15+080 obtained from our monitor- 
ing data in 2004-2006 dilfer from those determined by Schechter et 
al. (1997) and Barkana (1997) earlier. The differences for Tba and 
Tfic are well beyond the uncertainties reported in both publications 
and determined in the present work. While our time delay estimates 
for images A and C are rather close to the two previous ones, the 
delays for two other image pairs can not be regarded as consistent 
even marginally. 

• As could be expected, our estimates of time delays tac and 
Tbc result in larger Ho values than those reported by Schechter et al. 
(1997) with their estimates of time delays and with the ISXS, ISIS 
and ISIS+ models by Schechter et al. (1997). The new estimates of 
Hd are more consistent with the most recent Hq value obtained in 
the HST Key project (Freedman et al. 2001). 

• The new estimates of time delays in PGl 1 15+080 provide ad- 
ditional support for the family of models close to isothermal. As 
analyzed in details by Kochanek and Schechter (2004), the esti- 
mates of Hq with the use of time delay lenses are bounded by two 



limiting models: models with less dark matter (more centrally con- 
centrated mass profiles) produce higher values of than those 
with more dark matter. In particular, the constant mass-to-light ra- 
tio models set an upper limit on estimates of Ho, while the isother- 
mal mass distribution models are responsible for the lower limit of 
Hq. Our result is very important in this respect, since an isother- 
mal model is preferred for the lensing galaxy in PGl 115+080 for 
the reasons listed by Schechter (2005): (l)the velocity dispersions 
observed for an ensemble of lensing galaxies are consistent with 
the fundamental plane relations for ellipticals; (2) a majority of the 
nearby galaxies, as well as those lensing galaxies for which the ra- 
dial mass distributions can be measured, are very nearly isothermal. 
Since the PGl 1 15+080 lensing galaxy is by no means unusual, the 
isothermal hypothesis is most probable. 

In conclusion, recently published (Morgan et al. 2008) obser- 
vations of PG 1 15+080 in filter R during almost exactly the same 
time periods in 2004-2006 should be mentioned. We have used 
their table 3 photometry to compare to our light curves. The quasar 
brightness fluctuations which allowed us to determine the time de- 
lays are seen in their A1+A2 light curve quite well, but become 
undetectable in the B and C light curves because of a much larger 
scatter of the data points. 
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Table 6. Photometry of the PGl 1 15+030 images in filter V. 



Date 


JD 


Seeing (arcsec) 


Al 


A2 


B 




C 




20-04-2001 


2452019 


1.47 


17.025±0.097 


17.025±0.079 


18 


1.8±0.09 


IS 


i.232±0.088 


26-04-2001 


2452025 


1.258 


16.868±0.021 


17.19±0.017 


18 


1.736±0.02 


IS 


i.296±0.019 


27-04-2001 


2452026 


1.287 


16.874±0.037 


17.214±0.03 


18 


!.781 ±0.035 


18.298±0.034 


03-03-2002 


2452336 


1.171 


16.889±0.024 


17.252±0.019 


18 


!.782±0.022 


IS 


i.406±0.022 


08-03-2002 


2452341 


1.198 


16.985±0.046 


17.215±0.038 


18 


1.819±0.044 


IS 


!.435±0.043 


17-01-2004 


2453021 


1.038 


16.972±0.007 


17.185±0.006 


18 


1.782±0.007 


IS 


!.403±0.007 


13-02-2004 


2453048 


1.236 


17.0 ±0.02 


17.248±0.017 


18 


!.815±0.019 


IS 


i.424±0.018 


22-02-2004 


2453057 


1.262 


16.973±0.026 


17.28±0.022 


18 


!.843±0.024 


IS 


i.412±0.024 


26-02-2004 


2453061 


1.029 


17.04±0.014 


17.221±0.012 


18 


!.828±0.013 


IS 


i.418±0.013 


27-02-2004 


2453062 


0.845 


17.005±0.004 


17.243±0.003 


18 


1.831±0.004 


IS 


i.404±0.004 


28-02-2004 


2453063 


1.15 


17.021±0.019 


17.224±0.016 


18 


!.856±0.018 


IS 


i.41±0.018 


01-03-2004 


2453065 


1.059 


17.008±0.014 


17.237±0.012 


18 


!.832±0.013 


18.404±0.013 


10-04-2004 


2453105 


1.02 


16.934±0.006 


17.187±0.005 


18 


1.775±0.006 


IS 


i.353±0.006 


11-04-2004 


2453106 


1.183 


16.939±0.013 


17.187±0.011 


18 


1.792±0.012 


IS 


i.356±0.012 


14-04-2004 


2453109 


0.947 


16.947±0.005 


17.168±0.004 


18 


1.754±0.005 


IS 


i.352±0.005 


03-05-2004 


2453128 


1.203 


16.932±0.03 


17.196±0.024 


18 


!.699±0.028 


IS 


i.38±0.028 


09-05-2004 


2453134 


0.842 


16.96±0.004 


17.214±0.003 


18 


!.751 ±0.004 


IS 


i.405±0.003 


13-05-2004 


2453138 


0.842 


16.972±0.005 


17.218±0.004 


18 


1.772±0.005 


IS 


i.401±0.005 


17-05-2004 


2453142 


1.322 


16.942±0.026 


17.262±0.022 


18 


1.829±0.025 


IS 


i.437±0.024 


26-05-2004 


2453151 


1.169 


16.95±0.023 


17.279±0.019 


18 


!.792±0.022 


IS 


i.414±0.022 


04-06-2004 


2453160 


1.347 


17.077±0.027 


17.172±0.022 


18 


!.806±0.025 


IS 


i.438±0.025 


08-06-2004 


2453164 


1.233 


17.11±0.032 


17.145±0.027 


18 


1.845±0.03 


IS 


i.482±0.029 


05-01-2006 


2453740 


1.038 


17.316±0.008 


17.561±0.006 


19.09±0.007 


IS 


i.72±0.007 


08-03-2006 


2453802 


1.093 


17.313±0.011 


17.526±0.009 


19.083±0.01 


IS 


!.657±0.01 


15-04-2006 


2453840 


0.979 


17.2 ±0.01 


17.453±0.008 


18 


!.977±0.009 


18.576±0.009 



Table 7. Photometry of the PGl 1 15+030 images in filter I. 



Date JD Seeing (arcsec) Al A2 B C 



26-04-2001 


2452025 


1.197 


16.313±0.026 


16.698±0.021 


IS 


i.262±0.025 


17.781±0.024 


27-04-2001 


2452026 


1.144 


16.318±0.018 


16.699±0.015 


IS 


!.261±0.018 


17.781±0.017 


03-03-2002 


2452336 


1.01 


16.362±0.012 


16.654±0.01 


IS 


!.282±0.011 


17.817±0.011 


04-03-2002 


2452337 


1.351 


16.421±0.07 


16.578±0.058 


IS 


i.337±0.066 


17.845±0.065 


05-03-2002 


2452338 


0.858 


16.35 ±0.018 


16.662±0.014 


IS 


i.312±0.017 


17.793±0.016 


07-03-2002 


2452340 


1.288 


16.371±0.037 


16.695±0.03 


IS 


i.317±0.034 


17.841±0.034 


08-03-2002 


2452341 


1.139 


16.399±0.035 


16.661±0.028 


IS 


!.285±0.032 


17.836±0.032 


11-03-2002 


2452344 


1.043 


16.417±0.023 


16.66 ±0.019 


IS 


i.329±0.021 


17.853±0.021 


14-03-2002 


2452347 


0.928 


16.42 ±0.015 


16.67 ±0.012 


IS 


i.333±0.014 


17.828±0.014 


11-04-2004 


2453106 


1.136 


16.427±0.01 


16.622±0.008 


IS 


!.265±0.01 


17.814±0.01 


14-04-2004 


2453109 


0.773 


16.414±0.005 


16.623±0.004 


IS 


!.262±0.005 


17.801±0.005 


03-05-2004 


2453128 


1.125 


16.418±0.031 


16.598±0.026 


IS 


i.l82±0.029 


17.767±0.029 


09-05-2004 


2453134 


0.747 


16.427±0.005 


16.636±0.004 


IS 


i.25 ±0.005 


17.832±0.005 


13-05-2004 


2453138 


0.765 


16.428±0.005 


16.64 ±0.005 


IS 


!.252±0.005 


17.835±0.005 


17-05-2004 


2453142 


1.166 


16.426±0.021 


16.661±0.017 


IS 


!.265±0.02 


17.833±0.019 


26-05-2004 


2453151 


1.057 


16.437±0.016 


16.652±0.013 


IS 


!.253±0.015 


17.832±0.015 


04-06-2004 


2453160 


1.222 


16.464±0.054 


16.606±0.045 


IS 


i.213±0.051 


17.845±0.05 


08-06-2004 


2453164 


1.049 


16.526±0.027 


16.588±0.023 


IS 


!.256±0.026 


17.837±0.025 


05-01-2006 


2453740 


1.046 


16.672±0.009 


16.84 ±0.007 


IS 


!.509±0.008 


18.026±0.008 


08-03-2006 


2453802 


0.972 


16.675±0.008 


16.849±0.007 


IS 


i.512±0.008 


18.004±0.008 


15-04-2006 


2453840 


0.804 


16.61 1±0.006 


16.823±0.005 


IS 


i.413±0.005 


17.978±0.005 
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Table 8: Photometry of the PG 1 1 15+030 images in filter R. 



Date 



JD 



Seeing (arcsec) Al 



A2 



B 



20- 04-01 

26- 04-01 

27- 04-01 

03- 03-02 

04- 03-02 
08-03-02 
11-03-02 
17-01-04 

21- 01-04 

22- 01-04 
10-02-04 

13- 02-04 

16- 02-04 
22-02-04 

26- 02-04 

27- 02-04 

28- 02-04 
01-03-04 
30-03-04 

08- 04-04 

10- 04-04 

11- 04-04 

12- 04-04 

14- 04-04 

03- 05-04 

09- 05-04 

13- 05-04 

17- 05-04 

26- 05-04 

04- 06-04 
08-06-04 
08-02-05 

13- 02-05 

14- 02-05 

16- 02-05 

20- 02-05 

27- 02-05 

28- 02-05 

01- 03-05 

02- 03-05 

04- 03-05 

05- 03-05 

19- 03-05 

21- 03-05 

13- 04-05 

14- 04-05 

17- 04-05 

02- 05-05 

10- 05-05 
16-05-05 

20- 05-05 

03- 06-05 

05- 06-05 

06- 06-05 
20-06-05 

14- 12-05 

15- 12-05 
19-12-05 



2452019 
2452025 
2452026 
2452336 
2452337 
2452341 
2452344 
2453021 
2453025 
2453026 
2453045 
2453048 
2453051 
2453057 
2453061 
2453062 
2453063 
2453065 
2453094 
2453103 
2453105 
2453106 
2453107 
2453109 
2453128 
2453134 
2453138 
2453142 
2453151 
2453160 
2453164 
2453409 
2453414 
2453415 
2453417 
2453421 
2453428 
2453429 
2453430 
2453431 
2453433 
2453434 
2453448 
2453450 
2453473 
2453474 
2453477 
2453492 
2453500 
2453506 
2453510 
2453524 
2453526 
2453527 
2453541 
2453718 
2453719 
2453723 



1.226 
1.277 
1.126 
1.065 
1.359 
1.204 
1.266 
1.27 
1.45 
0.81 
1.254 
1.174 
1.451 
1.361 
1.065 
0.828 
1.225 
1.075 
1.157 
1.208 
1.033 
1.032 
1.145 
0.927 
1.134 
0.761 
0.827 
1.244 
1.005 
1.144 
1.138 
0.84 
1.192 
1.344 
1.461 
0.888 
1.011 
1.185 
1.232 
0.902 
1.017 
0.944 
1.097 
0.944 
1.168 
0.948 
1.016 
1.348 
1.161 
1.11 
1.208 
1.145 
1.381 
1.316 
1.421 
0.824 
1.111 
1.181 



16.625±0.011 

16.636±0.013 

16.63 +0.008 

16.693+0.014 

16.687+0.038 

16.691+0.049 

16.713+0.034 

16.762+0.031 

16.785+0.034 

16.767+0.004 

16.832+0.027 

16.799+0.019 

16.787+0.036 

16.832+0.037 

16.841+0.021 

16.799+0.004 

16.766+0.028 

16.813+0.018 

16.803+0.013 

16.772+0.01 

16.754+0.005 

16.735+0.006 

16.747+0.011 

16.754+0.004 

16.771+0.023 

16.78 +0.003 

16.781+0.004 

16.786+0.017 

16.781+0.013 

16.787+0.013 

16.862+0.023 

17.018+0.004 

17.06 ±0.011 

17.021+0.018 

17.08 +0.049 

17.042+0.007 

17.027+0.01 

17.038+0.015 

17.05 +0.013 
17.036+0.006 
17.044+0.006 
17.046+0.006 
17.057+0.008 
17.037+0.009 
17.017+0.01 
17.063+0.004 

17.06 ±0.011 
17.047+0.029 
17.101±0.012 
17.1 ±0.016 
17.101±0.027 
17.123±0.018 
17.225±0.041 
17.231±0.03 
17.188±0.104 
17.116±0.006 
17.12 ±0.008 
17.111±0.018 



17.069±0.009 

17.052±0.01 

17.071 ±0.007 

17.045±0.011 

17.079±0.031 

17.081 ±0.04 

17.081±0.028 

17.01 1±0.026 

16.998±0.029 

17.009±0.004 

17.015±0.023 

17.041±0.015 

17.087±0.03 

17.016±0.031 

17.007±0.018 

17.032±0.003 

17.072±0.023 

17.02 ±0.015 

17.013±0.011 

16.983±0.008 

16.985±0.004 

17.006±0.005 

16.991±0.009 

16.986±0.003 

16.983±0.019 

17.012±0.003 

17.02 ±0.003 

17.007±0.014 

17.041±0.011 

17.082±0.011 

17.013±0.019 

17.258±0.003 

17.26 ±0.009 

17.256±0.015 

17.268±0.041 

17.263±0.005 

17.322±0.008 

17.29 ±0.012 

17.258±0.011 

17.282±0.005 

17.272±0.005 

17.294±0.004 

17.289±0.007 

17.271 ±0.008 

17.324±0.008 

17.251 ±0.004 

17.306±0.009 

17.357±0.024 

17.323±0.01 

17.328±0.013 

17.367±0.022 

17.352±0.015 

17.257±0.035 

17.246±0.025 

17.329±0.087 

17.305±0.005 

17.294±0.007 

17.298±0.015 



18.573: 
18.579: 
18.579: 
18.59 + 
18.688: 

18.541: 

18.59 + 
18.574: 

18.612: 

18.596: 
18.62 + 

18.611: 

18.644: 
18.652: 
18.667: 
18.625: 
18.578: 
18.652: 
18.64 + 

18.615: 

18.597: 
18.598: 
18.578: 

18.591: 

18.567: 

18.601: 

18.598: 
18.579: 

18.601: 
18.615: 

18.667: 
18.83 + 
18.87 + 
18.842: 

18.941: 

18.868: 
18.867: 
18.89 + 
18.878: 
18.864: 
18.869: 
18.878: 
18.872: 
18.86 + 
18.852: 
18.848: 
18.872: 
18.866: 
18.892: 
18.904: 
18.943: 
18.944: 
18.965: 
18.964: 
18.985: 
18.922: 

18.915: 
18.911: 



=0.011 

=0.012 

=0.008 

0.013 

=0.036 

=0.046 

0.032 

=0.029 

=0.032 

=0.004 

0.026 

=0.018 

=0.034 

=0.035 

=0.02 

=0.004 

=0.026 

=0.017 

0.012 

=0.01 

=0.005 

=0.006 

=0.011 

=0.004 

=0.022 

=0.003 

=0.004 

=0.016 

t0.012 

=0.012 

=0.022 

0.004 

0.01 

=0.017 

=0.047 

=0.006 

=0.01 

0.014 

=0.012 

=0.005 

=0.005 

=0.005 

=0.007 

0.009 

=0.009 

=0.004 

=0.01 

=0.027 

=0.011 

=0.014 

=0.025 

=0.017 

=0.039 

=0.028 

=0.098 

=0.006 

=0.008 

=0.016 



18.095±0.01 

18.091±0.012 

18.106±0.008 

18.197±0.013 

18.215±0.035 

18.163±0.046 

18.178±0.031 

18.195±0.029 

18.186±0.032 

18.194±0.004 

18.205±0.025 

18.219±0.017 

18.206±0.034 

18.221±0.034 

18.216±0.02 

18.202±0.004 

18.177±0.026 

18.205±0.017 

18.16 ±0.012 

18.166±0.01 

18.163±0.005 

18.162±0.006 

18.148±0.01 

18.171±0.004 

18.169±0.021 

18.191±0.003 

18.201±0.003 

18.195±0.016 

18.218±0.012 

18.238±0.012 

18.255±0.022 

18.422±0.004 

18.442±0.01 

18.427±0.017 

18.439±0.046 

18.431±0.006 

18.438±0.01 

18.422±0.014 

18.433±0.012 

18.424±0.005 

18.419±0.005 

18.438±0.005 

18.435±0.007 

18.412±0.009 

18.431±0.009 

18.439±0.004 

18.448±0.01 

18.489±0.027 

18.492±0.011 

18.5 ±0.014 
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